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Introduction. The paper considers an axisymmetric problem of elasticity theory for a radially inhomogeneous 
transversally isotopic nonclosed sphere containing none of the 0 and zz poles. It is believed that the elastic moduli are 
linear functions of the radius of the sphere. It is assumed that the side surface of the sphere is fixed, and arbitrary 
stresses are given on the conic sections, leaving the sphere in equilibrium. The work objective is an asymptotic analysis 
of the problem of elasticity theory for a radially inhomogeneous transversally isotropic sphere of small thickness, and a 
study of a three-dimensional stress-strain state based on this analysis. 

Materials and Methods. The three-dimensional stress-strain state is investigated on the basis of the equations of 
elasticity theory by the method of homogeneous solutions and asymptotic analysis. 

Research Results. After the homogeneous boundary conditions set on the side surfaces of the sphere are met, a 
characteristic equation is obtained, and its roots are classified with respect to a small parameter characterizing the 
thickness of the sphere. The corresponding asymptotic solutions depending on the roots of the characteristic equation 
are constructed. It is shown that the solutions corresponding to a countable set of roots have the character of a boundary 
layer localized in conic slices. The branching of the roots generates new solutions that are characteristic only for a 
transversally isotropic radially inhomogeneous sphere. A weakly damping boundary layer solution appears, which can 
penetrate deep away from the conical sections and change the picture of the stress-strain state. 

Discussion and Conclusions. Based on the solutions constructed, it is possible to determine the applicability areas of 
existing applied theories and propose a new more refined applied theory for a radially inhomogeneous transversally 
isotropic spherical shell. 
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Introduction. One of the properties of materials affecting the stress-strain state of elastic bodies is their 
heterogeneity. Investigating the stress-strain state of inhomogeneous bodies based on three-dimensional equations of 
elasticity theory is associated with significant mathematical difficulties. 

A number of studies are devoted to the investigation of three-dimensional problems of elasticity theory for the 
sphere. 

In [1], based on the equations of elasticity theory for the sphere, a general solution satisfying the boundary 
conditions on the contour in the sense of Saint-Venant was obtained, the stress-strain state of the sphere was analyzed. 
In [2], based on the equations of elasticity theory for a thick isotropic sphere, homogeneous solutions depending on the 
roots of the transcendental equation are constructed. In [3], on the basis of solving three-dimensional problems of 
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elasticity theory for a sphere of small thickness, the accuracy of existing applied theories is studied, and a method for 
constructing refined applied theories is given. A three-dimensional asymptotic theory of a transversally isotropic 
spherical shell of small thickness is presented in [4]. An analysis of the three-dimensional stress-strain state of a three- 
layer sphere with a soft filler is described in [5]. In [6], the torsion problem is studied by the method of homogeneous 
solutions for a radially inhomogeneous transversally isotropic sphere of small thickness, when the elastic characteristics 
are changed by linear, quadratic and inversely quadratic laws along the radius. In [7], the torsion problem for a radially 
layered sphere with an arbitrary number of alternating hard and soft layers is studied. The existence of weakly damping 
boundary-layer solutions and a possible violation of the Saint-Venant principle in its classical formulation are shown. 
An applied theory of torsion of a radially layered sphere is constructed, which adequately takes into account the 
emerging features. In [8], the problem of elasticity theory for a radially inhomogeneous hollow ball is investigated using 
the finite element method and spline collocation. In [8], using the finite element method and spline collocation, the 
problem of elasticity theory for a radially inhomogeneous hollow ball is studied. The results obtained by finite element 
methods and spline collocation are compared. The axisymmetric problem of the theory of elasticity for a radially 
inhomogeneous transversally isotropic sphere of small thickness is studied by the method of asymptotic integration of 
the equations of the theory of elasticity in [9]. Inhomogeneous and homogeneous solutions are constructed. The nature 
of the stress-strain state is established. In [10], an axisymmetric problem of elasticity theory for a sphere of small 
thickness with variable elasticity moduli is considered by the method of homogeneous solutions. Asymptotic formulas 
for displacements and stresses are obtained, which provides calculating the three-dimensional stress-strain state of a 
radially inhomogeneous sphere. 

Materials and Methods. Deformation is considered within the framework of the linear theory of elasticity of a 
nonclosed sphere, whose material is transversely isotropic and inhomogeneous along the radial coordinate. The 
thickness of the hollow sphere is assumed to be small compared to the radius and size along the arc coordinate. 
Boundary conditions are considered that provide solving the problem in an axisymmetric formulation. We assume that 
the sphere does not contain any of the poles 0 and 7. In the spherical coordinate system, the area occupied by the 
sphere will be denoted by [= {r e[n;n], 9€[9;4,], d€[0; 2n]}. 

The linear dependence of the elastic properties of the material along the radius is considered: 

A, =a r, A, = ats A, =ayr, A), =ayr, A, =ayr, (1) 
where a\)),a\)’, a5)’, aS), a\) — some constants. 

The system of equilibrium equations in the absence of body forces in a spherical coordinate system 7r,0,¢ has 
the form [11]: 


do,, 100,, . 26, — Fy — Fog t O,,ctgO _ 








’ 0, 2 

or r 0@ r @) 

00,» : 16a ‘ 30,9 + (Fag — Tyg Clg _ 0, 3) 
or r oO r 


where, o 


rr? 


D, =U,(r,9), UV, =0,(r,9) as follows [4]: 


CF 45» — Stress tensor components, which are expressed in terms of displacement vector components 











o, =A, < i epee 0 4 ), (5) 

Oy = Ap a + 2] (4 + A,,)v, + A, v,ctg0 + A,, a (6) 
Og = Ary ” alae +A,)v, + Ayv,ctg0 + A,, a (7) 
oan ay(Sa-teo 8) " 


Substituting (5)-(8) into (2)-(3), taking into account (1), we obtain the equations of equilibrium in 
displacements. 
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Ou Ou 
by Bot +2¢ ae +2 (25). — bi) — bs ) Eu, + (20)? — bi! — bi) — bi? ) x 
Ou, Ou Oru Ou, Ou 
xe? | —2 4+ u,ctgd |+e(bO +b )| —2ctg6+ —4 |+ &6*b 2 4+—cte@ |=0, 9 
30 lS ) (bi: Sa dp & G00 “| ag? 36 & (9) 
o a Ou ou 
po) Me 9p 0 367, + e—* | 4.35) + +b Je? —2 + 
0p 0 1000 00 
é o Ou : , 
+e°bD oe ctg0 + - ) +ebo 5006 - (dS) ctg*O +b uy =0. (10) 
Here, p= nf =] — new dimensionless variable; ¢ = yn{ 2] a small parameter characterizing the thickness of 
E \N% ny 
the sphere: r, = . ia dine _Y, _Y% 0 4% ‘ . ge Hs 
phere; 7% = Ji mh; pel-Ll u, » Uy 3 b; G dimensionless quantities; G, some parameter 
1 1 0 


having the dimension of the elasticity modulus. 
We assume that the lateral part of the sphere boundary is fixed, i.e., 


=0, (11) 
up|... = 9. (12) 
We assume that at the ends of the sphere (on the conical sections), the stresses are specified 
Cwloa =f» Fpl. = Sr (P)- (13) 
Here, f,,(?), f5,(P) (s =1;2) — sufficiently smooth functions that satisfy the equilibrium conditions. 
Solutions (9), (10) are sought in the form [3, 4]: 
u,(P,9)=a(p)m(A); u,(p,4) = d(p)m'(A), (14) 


Where function m(@) satisfies the Legendre equation: 





m"(O)ctg0-m'(a)+{ 2 5 Jomo =0. (15) 
After substituting (14) into (9), (10), (11), (12), taking into account (15), we obtain: 


Ha p)+ Rea py +°| (40? 212 -269)-(2" +68 Jaco) 





(16) 
1 ' 
-[2-F}e[(o +99) 2%(o)~e(62 +08 +08) -299)a(p)]=9, 
" i, 1 
H9(a"(p)+2ed'p))-6"| 2-5) +(H9 09 +30) Jacoys in 
+6° (3b) +bi)) +B!) )a(p) + e(bL) +58) )a'(p) = 0, 
a(p) =0, upu p=. (18) 
d(p)=0, upu p=. (19) 
The solution to system (16), (17) has the form: 
a(p)=e” [ pe’C, + pe’C, + pe", + pe Pd,| > (20) 
d(p)=e” [te Gage "C462" C, +96 YC, Al : (21) 


where C, (n= 1,4) — arbitrary constants, 


1 
Py = Oia Sy [2 -7 |e +(b!) — bY) — 4b? ) ; 





= (0), (0) (0) 4. p(0) 1 7,00) (0) \. 
t, =— (bf) +b) )s, -(2b0) +b) +8) -BY)); 
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C= (D\ bY |S; (2b) +b +b +b —b? : Ss, — equation roots 
1 
by Os + IG = ra: >? + 26D BD) -B0BO )+ BOY — SBD EO + 4B0 DY - 
= 26.5.9 — 52 268]? + IG 7 a ~ 4b +25 +26 + : (22) 


«(2 -F ag vag ag ane |-2- + (202 +08 +00 - -1®) =. 


The system of linear algebraic equations with respect to C,,C,,C,,C,, is obtained through satisfying the 


homogeneous boundary conditions (18), (19). The zero equality of the determinant of this system is a condition for the 
existence of nonzero solutions and leads to a characteristic equation with respect to spectral parameter 


A(z;é) = (p,q — Pr \(4P2 —4P, )sh? (e(s, +5, ))+ 
+ (Dit, — Prd )( P14 — Pot, ) sh” (e(s, —s,)) = 0. 


Equation (23) has countable set of roots Z,. The general solution to the problem is obtained by summing over 


(23) 


the roots of equation (23) 


u, = >\M,a,(p)m,(0), (24) 
uy = Md, (pm! (0), (25) 


where 
a,(p)=e” [21 (ea, —e gp, )+ P2 (e°*e, =e OO, )|. 
d.(p)=e” [nea —qe’a, +t,e°" a, — ne*”a, |, 


€(2s)-8)) 


a, = pia (> a t, (DD Os (pt, — Pr, Je “os 





él (5, +255) 


@, = ppt, (g,-t, )e™ -1, al pot, eo") +, (pity — pot, Je 


ae + PQ (4, —f, )e* > 


=a (Pit, — Pon )e ~q (pt, —P2 t, )e* val + pf, (9, -14,)e™. 


The set of roots of equation (23) at ¢ +0 consists of countable sets of roots 


@, =t, (Dd, — Dig Je" — 4, (D,da — Pat, Je 


&(8,—-25;) (25, 





= Oe 4 Ole). (26) 
é 
For 6),, we have: 
1°. At b, >0, b?-b, >0: 
it —s,)(b? —b s,s, )sin((s, +58,)5)+(s,+8, )(ay? +bs,s, )sin((s, —s,)6)=0, (27) 
where 
5, =fb +b? =; 5, = Jo - Jb? -4,; 
b, = (20H) ne +OOY -BP BD); B, =BOaPY. 
2°. At b >0, b?—b, <0: 
| 20(.a? - (bi) (a? - B°)—b4)) | sh(26er) + 
(28) 
+a 26) 6? + (Bi? (a? - 6°) a.) ) |sin(208) = 
where 
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3 Ath >0, & =b,: 


(Bis? —bf )sin(2sd) +2( b's” +b) )ds =0, (29) 
where s = lB 
4°. At b, <0, b?—b, >0: 
(3,3 Non —bs,s, )sh((s, +5s,)5)+(s, +8, )(a? +b s,s, )sh ((s, -s,)5) =0, (30) 


where 
s, = yb |- Je -4,; 5, = yfo|+ Jb? -5,. 
5°. At b, <0, b—b, <0: 
B| 20a? -(b\) (a? 6°) -_b.) ) |sin (26er) + 


(31) 
tal 20) 6° + (H{ (a? —B°)-Bg)) | sn (268) =0, 
where 
a= Vio |-V? -2, =1(@-if); 5, = lo l+Vb?-2, = +(at+if). 
6°. At b, <0, bf =b,: 
(2's? — bi )sh(2sd) + 2(d(s° +b )os =0, (32) 


where s =,/|b,|. 


Equation (27)—(32) has a countable set of solutions. 
Let us present an asymptotic construction of solutions corresponding to different groups of roots of the 
characteristic equation (23). Substituting (26) into (24), (25) and expanding the resulting expressions in powers €, we 


have: 


le 
a) u, )= des {(o2 +b)| (os) +b”s; Js cos(6,s, )sin(d,s,2)— 
k=1 
-(b? +505? is cos(6,s,)sin(d, sp) |+O(e)\m, (A), (33) 
uy(P = YE {(a2 + BOs: (os +()s; )[ cos( (5,5, )cos(6,s,2)—cos(5,s, )cos(6,8,p) ]+O(é e)}m 1 (9), (34) 


k=l 


where 6,, are the solutions to equation 


(5, —s,)(b? pe +58,)5)+(s, +5,)(b +b())s,s, )sin((s, —s,)5)=0 (35) 
6) u, )= BP 6: {(o? +b.))[ (oe? +b\)s; )s, sin(6,s, )cos(d,s,p) - 
~(b) +4,s?)s, sin(6,s, )cos(5,s,) |+O(«)}m, (8), (36) 
uy (p )= Vest {(o? +b\s; (by) +b\)”s; )[ sin( d,s, )sin(d,s,) —sin(d,s, )sin(d,5,p) |+O(e) (2)}m (8), (37) 


k=1 


where 6,, are the solutions to equation 


(s,—s, )(ay? —bs,s, )sin((s, +8,)5)-(s, +5, y(ay +b s,s, )sin((s, —s,)5)=0. (38) 


99 
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a) u,(p;0)= > £53 (dy) +b) ){[ Bsin(5, Be) sh(d,ap) — acos(6,8p)ch(5,ap) |: 
| 2086) cos(6,)sh(6,a)+ (a? (a? —£ ) —b )sin(6,8)ch(5,c) | - 


- [a sin(d,8)sh(6,ap)+ Bcos(6,fp)ch (d,ap) | . 


| 2ap0\? sin (5,B)ch(5,a)—(b\? (a? - 6°) —b.) cos (5,8) sh (5,2) + o.e)| m,(0), (39) 


u, (p39) = > ES! {[ 200 sin(6,8e)ch(6,a@p)— (a? (2? —£ ) —b\ ) “ 
xeos( 5,2) sh(6,ap) ]| 2aBb\) cos(6,8)sh( 6,02) +(b\) (a? - 6?) —bi? )sin (5,8) ch(,cr) |= 
-| 2b cos(d8,p) sh(6,ap)+ (a) (a? —£ ) —b )sin (6,80) ch (6,2) | x 


| Zapp)? sin (6,8) ch(6,2)—(8\? (a? — B°)—b,? )cos(6, B)sh( 5,2) | +0(«)} m, (8), (40) 


where 6,, are the solutions to equation 


Bl 2)er — (6 (a? - B?)-b ) [sn (250) + ar| 26? B +(? (a? - B°)-Be )] sin(258)=0. (41) 
6) u, (2:0) = SE 5 (2: +. ){[ Bsin(5, Bp) ch(6,ap) — arcos( 5, Bp) sh( 5,ap) |: 
| 206? cos(5,8)ch(5,a) +(B (a? - B*) bs )sin(6, B)sh(6,c) | 7 
—[ asin(5,8p)ch(6,ap) + Bcos(d, Be) sh(6,ap) || 2aBb{)” sin(d, 8) sh(d,a) - 
-(b\? (a? ~ B°)B8) Jeos(5,8) ch(5,a)]+O(2)}m, (8), (42) 
u, (p39) = > Es {[ 22800 sin(6, fp) sh(6,a@p)— (a, (a? —£ ) —b )x 
xeos(5,2)ch(5,ap) | 2aBb\ cos( 6, 8)ch( 5,01) +(b\? (a? — 6?) —b4) )sin(5,8)sh(d,a) |- 
-| 2b) cos (dB, p)ch(6,a@p)+ (ay (a -£ )-o? )sin(5, 8p) sh (6,2) | x 
x| 2a,p0\? sin (6,8) sh(6,a)—(b\) (a? — B°)—b.) )eos (5,8) ch(5,c) +0(«)} m\ (0), (43) 
where 6,, are the solutions to equation 
B| 2b ee — (5° (a - B?)-b )] sh(26e) -a [ 25° B +(B (a? -B°)-10 ) |sin(208) =0. (44) 
a 


a) u,(p;0)= > {(a0s? +b\!))(cos(5,s)cos(d,sp) + psin(d,s)sin(d,sp)) + 
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(ays? — be) 
+ ‘sin (d,s )eos(6,sp) + O(¢) m, (0), (45) 
SO; 
in 0) (0 
0)= ys" (bis? +b.) x 
il (B{))s? +b!) 5,s 
x(psin(d,s)cos(d,sp)—cos(6,s)sin(d,sp)) +O(e)} m. (0), (46) 
where 6,, are the solutions to equation 
(B's? bf )sin(2sd) +2( B's? + bf) )ds = 0 (47) 


2 BOs? — pO 
6) u,(p;0) => E,” eo (5,5) sin(5,sp)-(6)'s? +b\))x 


k=1 iS 


x(sin (6,5) sin(d,sp) + p.cos(6,s)cos(d,sp))+ O(e)} m, (0) (48) 


. BOs? +B) 
‘g\= Eo ( i 44 
uo(7:9) ie (6 +5) 5,5 


x| pcos(6,s)sin(5,sp) —sin(d,s)cos(6,sp) |+ O(e)} m,(@), (49) 


where 6,, are the solutions to equation 
(2)? —bit) )sin(2sd) -2(d\)'s* +i) ds =0 (50) 
4°. In case b, <0, b> —b, >0, asymptotic formulas for displacements are obtained from (33)—(38) through 
replacing s,,s, with is,,is,. 
5°. In case b, <0, b —b, <0 , asymptotic formulas for displacements are obtained from (39)—(44) through 
replacing s,,s, with is,,is,. 
6°. In case b, <0, b> =b,, all asymptotic formulas for displacements are obtained from (45)-(50) through 


replacing s withis . 
For roots (26), the main term of the asymptotic solution of equation (15) at ¢ +0 takes the form [9, 10]: 





Eel 8% (8 4) |(1+0l6)); near 0=6,, 











m,(@) = (51) 
exp| ‘J -6,, (@-6 | 1+O(e€)); near 0=86,. 
a wea 0k ( , ) ( ) 2 
We represent displacements in the form: 
u,(p,) =) Ea (p)m, (9), (52) 
k=l 
Ug(p,0) = ¥ E,d,(p)m, (8). (53) 
k=1 
We represent stresses 0, and o,, in the form: 
00 = DE, (a2 (p)m, (9) + of? (pm, (O)ctg), (54) 
k=1 
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O59 = DEO (P)n, (9), (55) 
kal 
here, 
(1) 1 (0) (0) (0) (0) 2 1 : 
O71 (P) = = [bs a,(p)+ e (BS +b;; Ja, (p)—ébyy’| Zz — Z d,(p) |; 


of? (p) = (BS) —b )d, (p): 
(0) 


ox(p)=-[di(p) + 6(a,(0)-d,(p)} 


The character of solutions (33)—-(50) depends essentially on the type of roots 6,,. The first boundary-layer 
terms of these solutions correspond to the Saint-Venant edge effect [4]. In the case of imaginary roots 6,,, these 
boundary layers have weak damping. Thus, the stress-strain state is far enough away from the ends and significantly 
depends on them. That is, in this case, the transversely isotropic properties of the inhomogeneous material significantly, 
in comparison to the isotropic material of the sphere, change the pattern of the stress-strain state. At the same time, for 
real or complex 6,,, the pattern of the stress-strain state of an inhomogeneous sphere for such materials qualitatively 
coincides, differing in the decay rate of the above-described Saint-Venant boundary layer solutions of an 
inhomogeneous plate. 

From (51), it turns out that when moving away from the conic sections 6=6,(j =1,2), solutions(33)(50) 
decrease exponentially. 


Since the constructed solutions satisfy the equilibrium equation and boundary conditions on the side surface, 


the Lagrange variational principle takes the following form [4, 11]: 





> [[les =Ji (p)) du, + (O45 — f,;(p)) du, | bs edo =. (56) 


j=l 
Substituting (52)-(55) into (56) and taking dE , aS independent variations, we obtain an infinite system of linear 


algebraic equations 
Ede =F, G=1L2,...), (57) 
k=l 

Here, 


1 2 1 
Gn = | oP(p)4, (nied Som, (0, )m’, (| + [ of (p)d,(p)e*?d px 
=] s=l =I 


«(Som (8,)m(8, v0, | fox(o)a, (reap Som (,)m, @)), 


2 


a=) m 0.) f.(p)d,(p)e”dp +m,(0,){ f,,(p)a, (pea | 


s=l 

System (57) is always solvable under physically meaningful conditions imposed on the right side (57). The 
solvability and convergence of the reduction method for (57) is proved in [12]. 

Using the smallness of parameter €, it is possible to construct asymptotic solutions of system (57). 

Research Results. The structure of the stress-strain state of a radially inhomogeneous transversally isotropic 
sphere of small thickness is analyzed under kinematic conditions on the side surface. It is shown that, in the case of 
fixing the side surface, the character of the solution is determined by the boundary layers. It is found that the asymptotic 
decomposition of the stress state starts with a solution describing the Saint-Venant edge effect in the theory of 
transversally isotropic inhomogeneous plates. In the case of transversal isotropy of the radially inhomogeneous material 
of the sphere, some boundary layer solutions decay very weakly, they can penetrate deep far from the conic sections and 
change the pattern of the stress-strain state. Asymptotic relations for displacements and stresses are derived, which 


provide calculating the three-dimensional stress-strain state of a radially inhomogeneous transversally isotropic sphere 
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of small thickness with any predetermined accuracy. It is shown that the root branching generates a countable set of 
new solutions for a transversely isotropic radially inhomogeneous sphere. 

Discussion and Conclusions. An asymptotic analysis of the stress-strain state of inhomogeneous shells, based 
on three-dimensional equations of elasticity theory, makes it possible to establish the limits of application of 
approximate theories. The identified behavior pattern of the solution far from the ends for different boundary conditions 
on the side surfaces can become the basis for creating refined applied theories for calculating the deformation of a 
radially inhomogeneous transversally isotropic spherical shell of small thickness. One of the applications of the 
asymptotic analysis performed can be the calculation of shells with thin coatings, in which, in this case, a radial 


inhomogeneity arises [13, 14]. 
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